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Bose-Einstein Condensation with Entangled Order Parameter 
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We propose a practically accessible non-mean-field ground state of Bose-Einstein condensation 
(BEC), which occurs in an interspecies two-particle entangled state, and is thus described by an 
entangled order parameter. A suitably defined entanglement entropy is used as the characterization 
of the non-mean-field nature, and is found to persist in a wide parameter regime. The interspecies 
entanglement leads to novel interference terms in the dynamical equations governing the single 
particle orbital wavefunctions. Experimental feasibility and several methods of probe are discussed. 
We urge the study of multi-channel scattering between different species of atoms. 
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For a mixture of two distinct quantum many-particle 
systems, the usual approach is to associate an order 
parameter to each. For example, in a two-component 

BEC USUI!, 

each component is separately described 
by a condensate wavefunction, although the two wave- 
functions are classically coupled. Beyond this approach 
is BEC of spin-1 H 0, 0, @ or spin-± atoms [1 GjJ, 
but under consideration is one many-particle system in 
which all the atoms are indistinguishable and each atom 
can jump in between different spin states. In this Letter, 
we propose a peculiar mixture of two distinct species of 
atoms a and b, which represents a novel situation different 
from both cases above. Each atom can flip between two 
pseudo spin states f and j, but cannot transit between 
a and b. The total occupation number of each species 
is conserved, but that of each spin state of each species 
is not conserved because of scattering between different 
species. It is shown that the ground state is BEC of 
inter-species pairs in a two-particle state which, as the 
global BEC order parameter, is quantum-entangled. As 
such, this work extends the notion of entanglement to 
order parameters of many-body systems, and represents 
a direction different from other works on entanglement 
in BEC I n particular, the entanglement between 

the two species can serve as a characterization of the 
non-mean-field nature. The detection methods are dis- 
cussed. We also derive the dynamical equations of the 
single particle orbital wavefunctions, each containing an 
interference term which is absent in BEC mixtures stud- 
ied previously. As a condensate of nonlocal Bose pairs, 
our system is also related to the condensation of bound 
molecules and of paired fermions [l2T |. Experimentally, 
the two pseudo spin states can be implemented as two 
hyperfme mj? states in a magnetic trap. The atoms can 
be constrained in the two pseudo spin states by exploit- 
ing the conservations of energy and of the total mp in 
a scattering. Because of the implemental feasibility as 
well as the persistence of the entanglement in a wide pa- 
rameter regime, this BEC system represents a practically 



accessible non-mean-field ground state, in contrast with 
previous BEC models. 

In ignorance of the depletion, the orbit of each 
atom, of species i (i = a, b) and with pseudo 
spin er, is constrained in the single-particle ground 
state <j)i a {Yi). The Hamiltonian can be writ- 

ten as H = £ CT f itr N ia + § Y. aa > K 3 N f N ™' + 
T,**'Klf)N a(r N b(T , + ^( t 0i 6[6 T + a\a^b\b{), where 



a a and b a are annihilation operators, N aa 
N b<T — b\b a , K^/'s are related to scattering lengths [l 
fi a = ti a — K„a /2, where ti a is the single particle en- 
ergy. e aT - tai = fb[ - e&t- N i = N 'i.1 + N n is a constant. 
The last term in TL describes the spin-exchange scatter- 
ing, which changes the occupation of each spin state of 
each species and causes interspecies entanglement. De- 
fine S a = Y^a,a> a i S <T<r' a o-' and S 6 = Ha,a' b t s <T<T' b <7' , 

where s is the spin-^ operator. The common eigenstate of 
S% and S az is \S a ,m a ) = +m a ) a] \^ -m a ) ai . Thus 
we transform the Hamiltonian into that of two coupled 
big spins S a = N a /2 and S b = N b /2, 
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where J z 



K, 



~jf(SaxS bx + S a yS b y) + S aZ S bz 

+B a S az + B b S bz + C a S az + C b S bz 
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- K^ b) )/2, B b J z = - f bl + N b (K\ bb) 
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(66), 



U - "Ti - -"-J.T )/ 2 - E ° = + 
(1/8)E,(E^4S)^ + (m(T,^Kif))N a N b is a 
constant and is neglected henceforth. 

Ti. conserves the total z-component spin S z — (JV -f — 
N a { + N b] - N H )/2. The total spin is S = S a - 
S b ,---,S a + S b , letting N a > N b . For a given 
S z , any eigenstate of Ti can be written as |^) = 
g(m)\S a ,rn)\S b , S z — m), with the summation range 
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max(—S a , S z — Sb) <m< min(S a , S z + Sb)- 

At the isotropic point, i.e. if K e = J z while B a = 
Bb = C a = Cb = 0, the Hamiltonian is simplified 
as H = J 2 S a • Sb, whose eigenstates are nothing but 
IS*, S z ), for which g(m) is the Clebsch- Gordon coefficient 
(S a , m, Sb, S z — m\S, S z ). The eigenvalues are degenerate 
for different values of S z . There are N a — Nb + 1 degener- 
ate ground states \Gg K ) = \S a — Sb, S z ), which are found 
to be exactly 

\G s ,)=A(a\r(a\r±(a\b\-a\b\) N >\0), (1) 

where n T = N a /2 - N b /2 + S z , n l = N a /2 - N b /2 - S z , 
A = [(N a - N b + l)\/(N a + l)\N b \(N a /2 - N b /2 + 
S z )\{N a /2 - N b /2 - S z )\} 1/2 is the normalization fac- 
tor. Particularly, for N a = Nb = N, the ground state 
is uniquely |G ) = (■^/W+1N\)- I (a\b\ - a\b\) N \0). 

The non-mean-field nature of these ground 
states can be well characterized by the entangle- 
ment between the two species. Rewrite \Gq) as 
( % ^VTT)- 1 ELo(-l)^"V)aTl^-^)al|iV-m) fcT |m) bi 
The entanglement between the occupation of each one of 
the four single particle states and the rest of the system 
can be quantified as the von Neumann entropy of its 
reduced density matrix 0], which is log Ar+1 (iV+ 1) = 1, 
where we set the logarithmic base to be TV + 1 because for 
each single particle basis state, the occupation number is 
(N + l)-valued. In a sense \Go) is maximally entangled, 
due to equal superposition of the N + 1 quart-orthogonal 
states. As each species comprises two spin states, it can 
be seen that the entropy of the entanglement between 
the total a species and the total b species is given by 
the same entropy. For the general eigenstate \ty) above, 
one can also write it in terms of occupation numbers in 
four single particle states, and obtain the entanglement 
entropy as — ^ m g 2 {m)\og Nb+1 g 2 (rn), which must be 
less than 1 if \^) ^ |Go), as g(m) 2, s are not all equal for 
different m's. 

\Gs z ) is a condensation of interspecies 
pairs in the same two-particle entangled state 
</>(r Q ,r ) = ^[<?Mra)| t)a<hl(n)\ l)b - 
4>ai(r a )\ 4>a0bT( r fc)l T)b], plus the extra a-atoms. 
This can be seen by noting (ajfrj — ajoj)^ 6 = 

[V2 J d 3 r a d 3 r b ^i{Ya)'ilji{Yb)<t>{Ya,n)\ Nb , where 

VvM = Y,a a <r ( t ) aa{ra)W)a, ^b( r ) = J2a b cr(/>ba( r b)\^)b- 

Thus the BEC order parameter, given by </>(r Q ,r b ), 
is quantum-entangled [lflj . Like Cooper pairing, 
the formation of the entangled two-particle state 
lowers the energy, as can be seen by consider- 
ing one a-atom and one &-atom, with Hamiltonian 
fe(ra)+/»(r6)+^i(r a -r 6 )+^ 2 (r -r 6 )(| U)(U l + l ITXlt 
|), where C/ 2 > 0. Compared with the separable state 
4>a{ r a)4'b( r b)\&)\< jl ), the entangled two-particle state 
)(j>b ( r b)( I T)| I) — I 1)1 ]))/V2 lowers the energy 
b y I U^ ab \r a - r b )|(/> Q (r a )0 b (r 6 )| 2 . For such energetic 
reason, our BEC occurs globally in the interspecies 



entangled pair state </>(r a ,r b ), rather than separately in 
the two species. 

This entanglement can be experimentally probed 
by several means. First, because of entanglement, 
the particle number of each spin state of each 
species is subject to strong fluctuation. For \Gq), 
(Nur) = N/2, (N 2 a ) - ( N la ) 2 = N(N + 2)/12, 
and thus ^/(N 2 a ) - (N ar7 ) 2 / (N a<7 ) « l/y/3. Note 
that the quantum mechanical average can only be 
obtained through several runs of the measurement. 
Nevertheless, as the density is a self-averaging quan- 
tity |l6| . one can study, in a single image, the 
fluctuation of the density pj a fa) N ia \(j)i a (ri)\ 2 . 
In terestingly, ^ {p ilJ (Y i ) 2 ) - (p^fa)) 2 /^^)) = 

V(ND ~ <^)7<AW- 

Secondly, this entanglement is also indicated by 
nonvanishing of the connected correlations between 
particle numbers of the two species, G a ^ a i = 
{N aa N h a') - (N aa )(Nba')- For \G ), C a , a = 
-N{N + 2)/12, C a ,a = N(N + 2)/12, where 
a 5^ a. Again, these connected correlations can 
be measured through the corresponding quantities of 
the densities, g(r a , a; r b , a') = (/wfaW' (r b )) - 

(Paa(r a ))(pba'(rb)) = C a ^\<\> )(f>ba(^b)\ , and 

5( r a, v\ r b , a')/(p aa {r a )) (p ba > (r b )) = C ai<r , / (N aa ) (N ba >). 

These density fluctuations and correlations are those in 
the original condensate, so it is not necessary to switch 
off the trap to let the condensate freely expand. Never- 
theless, one can do so in order to obtain a larger image, 
which equally contains information about the entangle- 
ment. This is because the free propagator does not af- 
fect the entanglement. Moreover, the density fluctuations 
and correlations in the expanded gas are proportional to 
corresponding quantities in momentum space [lr3 | , which 
are related to those of the particle numbers in a way sim- 
ilar to that in the coordinate space, as described above. 

Thirdly, this entanglement can also be probed by 
detecting atoms which leave the trap. Consider 
two-atom measurements in which the joint probabil- 
ity is obtained. In measuring the spin state of 
an a-atom, the probability of obtaining a is Pj CT = 
(a\a a ) / J2a< ( a l ,a o-') [HI- For a joint measurement of 
the spins of an a-atom and a 6-atom, the probability 
of obtaining a for the a-atom while a' for the 6-atom 

is P a , a i = (bi,ala a b a ,)/J2a a ,<? b ( b l b a t a a <>a b vb)- A non- 
entangled, or mean-field type, state is in the form of 

, , * Ni i N 2 4- N 3 + Na, , 

(ViV 1 !iV 2 !iV3!A4!)- 1 at i aL fl &t bl & |0), where n 
and m are two arbitrary "directions" . Hence for arbi- 
trary chosen a a and <r b , Pa a ,<r b — P<j a Pcr b - In contrast, 
for an entangled BEC considered here, P Ga ,a b ^ Pa a P<r b - 
Indeed, for |G ), P„ a = P ab = 1/2, but P n = P iT = 
(2N + l)/6N while P TT = P u = (JV - 1)/6N. 

Finally, this entanglement has feedback effects on 
the single particle orbits. We can derive the equa- 
tions governing the orbital wavefunctions, as the coun- 
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FIG. 1: The entanglement as a function of K e /J z , with 
B a — Bb — C a = Cb — 0, S z — 0. The entanglement measure 
used both here and in the next three figures is the von Neu- 
mann entropy of the reduced density matrix of either of the 
two species, calculated by considering the occupation entan- 
glement [14j . 

terparts of the Gross-Pitaevskii equation, using the 
method of Ashhab and Leggett [lfj. The minimiza- 
tion of the expectation value of TL in the ground 
state \G ) gives rise to {-^V 2 + £7 aCT (r) + [JV(2V - 

l)/3]g { *l a) \<t>a*(r)\ 2 + [N(N- l)/6]g\f \$ aS (r)| 2 + [N{N— 

l)/6] ff ^ b V^(r)| 2 + [iV(2iV+l)/6].g^ ) |^(r)| 2 }0 Q(T (r)- 
[N(N+2)/l2]g e 4>l-(r)<t> ba (r)4> a , 9 {r) = fj. aa (t> a a(r) and the 

Hi) 

likes for the other three wavefunctions, where q„J and 

' J (7(7 

g e are the coefficients in front of the integrals in K^) 
and K e |l3l |. The term proportional to g e is an inter- 
ference term, which persists even in the Thomas-Fermi 
limit. It is absent in previous models and may lead to 
novel physical properties. 

We now study how this entanglement survives the cou- 
pling anisotropy J z ^ K e and the nonvanishing of B a , 
Bf,, C a and Cb- We use the Lanczos method to numeri- 
cally determine the ground states for different parameter 
values. First we study how the entanglement varies with 
K e /J z , setting B a = Bb = C a = Cb = 0. As shown 
in Figs. 1 and 2, the entanglement is maximized when 
K e = J z , and is maximized at S z — 0. This maximum is 
1 if S a = Sb, and is less than 1 otherwise. Roughly speak- 
ing, the larger the spins, the sharper the deviation of the 
entanglement from the maximum, but the details depend 
on all of S a , Sb and S z . The decrease for K e /J z > 1 is 
much slower than the increase for K e /J z < 1. 

Under typical values S a = 12000, S b = 10000, S z = 
1000, K e / J z = 1.2, we study how the entanglement varies 
with C a and Cb, while B a = Bb = 0. Because J z , C a J z 
and CbJ z are about of the same order of magnitude of 

(ii) 

Kp£f s, we consider the range —2 < C a < 2 and —2 < 
Cb < 2. As shown in Fig. 3, the change of entanglement 
from the case of C a — Cb = 0, in which the entanglement 
is 0.6, is quite limited. 

Then we typically choose C a = 0.2 and Cb = 0.4, and 
study the variation of the entanglement with B a and Bb- 
As B a and Bb can be as large as the order of magni- 



FIG. 2: The entanglement as a function of K e /J z , with B a = 
B b = C a = C b = 0, S a = 12000, S b = 10000. 




FIG. 3: The entanglement as a function of C a and Cb, with 
S a = 12000, S b = 10000, S z = 1000, Ke/J z = 1.2, B a = 
B b =0. 

tude of N a and Nt, or even larger, we consider the range 
-10000 < B a < 10000 and -10000 < B b < 10000. Fig. 4 
indicates that in such a parameter range, the variation 
of entanglement, with the maximum 0.57, is still quite 
limited. 

These results show that BEC with an inter-species 
entangled order parameter, as a highly non-mean-field 
ground state, can occur in a quite large parameter 
regime. The reason why the entanglement is not eas- 
ily destroyed by anisotropy and C a , Cb, B a and Bb is 
that the entanglement is caused by the exchange inter- 
action. This leads to more relaxed condition for the per- 
sistence of entanglement than for the energetic pertur- 




FIG. 4: The entanglement as a function of B a and Bb, with 
Sa = 12000, Sb = 10000, S z = 1000, Ke/Jz = 1.2, C a = 0.2, 
C b = 0.4. 
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bation. Moreover, around C a « — C b , there is a can- 
celing effect in the additional term in the Hamiltonian, 
C a Saz + @bSbz = + Cb)S az — 2C b S z S bz + C b S z . Sim- 
ilarly, around B a rs £?;,, there is a canceling effect in 
5 a S az + B b S bz = {B a - B b )S az + B bz S z . The cancel- 
ing effect further reduces the change of entanglement, as 
exhibited in the Figs. 3 and 4. 

It is important that the non-mean-field ground state, 
with inter-species entanglement, can be accessed in ex- 
periments. A crucial issue in experimentally accessing 
the non-mean-field ground state is the energy gap. At 
the isotropic point of our Hamiltonian, in analogy with 
BEC of spin-1 bosons, the gap is of the order of ~ 1/V, 
thus the ground state is very difficult to be reached in 
practice |6j. This is spontaneous symmetry breaking. 
However, it can be estimated that when B a — B b is of 
the order of ±Y or C a + C b is of the order of —1, the 
energy gap is of the order of N/V, which is finite in the 
thermodynamic limit. On the other hand, we have seen 
above that in this parameter regime, as far as K e is larger 
than or not much smaller than J z , the entanglement is 
significantly nonvanishing. Even if B a , B b , C a and C b are 
not so large, the finite energy gap can still be opened up 
when J z 3> K e . But then the entanglement tends to van- 
ish. One may resort to adiabatic switching, which might 
be realized by using Feshbach resonance. One starts with 
the gapped case J z ^> K e . By adiabatically tuning K e 
towards K e w J z , the ground state evolves into a highly 
entangled one, though with vanishing gap. 

The hyperfine-Zeeman energy levels depend only on F 
and rriF, not on the species [8j. Therefore we can simply 
represent the same pseudo spin state of the two species 
by the same actual hyperfine state. To exactly realize our 
model, only forward and exchange scattering channels are 
needed, thus the scattering into other hyperfine states 
needs to be suppressed. A choice is | f) = \F = 2, nip = 
2); I I) = \F = 2,m F = 1). One may also use | |) = 
\F = 2,m F = 2), | |) = \F = l,m F = 1), or use | T) = 
\F = 2,m F = -1), | i) = \F = l,m F = -1), which are 
the spin states in the well-known two-component 87 Rb 
BECs 0. In these previous experiments, because the 
particle number in each spin state is conserved, the BECs 
are equivalent to a mixture of two species, each in a fixed 
spin state, which was implemented in 41 K and 87 Rb, both 
in \F = 2,m F = 2) jj]. These experimental experiences 
can exploited to realize our model. For the latter two 
choices of spin states, as the energy difference between 
the two spin states with different F's is very large, B a 
and B b can be of the order of N or larger. Likely this can 
produce the gap. Lack of data on interspecies scattering 
prohibits more detailed realistic discussions here. 

To summarize, we have proposed a practically accessi- 
ble non-mean-field BEC, which is a condensation of inter- 
species spin-entangled pairs and can be characterized by 
a suitably defined entanglement entropy, which is max- 
imized in the ground state at the isotropic limit, but 



remains significant in a wide parameter regime. We have 
also derived the dynamical equations of the single particle 
orbital wavefunctions, where there are novel interference 
terms. This non-mean-field ground state can even be 
gapped. We have discussed the experimental feasibility 
and put forward several means of experimental detection. 
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